Appearance of the central singularity in spherical collapse 
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Abstract 



We analyze here the structure of non-radial nonspacelike geodesies terminating in 
the past at a naked singularity formed as the end state of inhomogeneous dust collapse. 
The spectrum of outgoing nonspacelike geodesies is examined analytically. The local 
and global visibility of the singularity is also examined by integrating numerically 
the null geodesies equations. The possible implications of existence of such families 
towards the appearance of the star in late stages of gravitational collapse are considered. 
It is seen that the outgoing non-radial geodesies give an appearance to the naked 
central singularity as that of an expanding ball whose radius reaches a maximum 
before the star goes within its apparent horizon. The radiated energy (along the 
null geodesies) however, is shown to decay very sharply in the neighborhood of the 
singularity. Thus the total energy escaping via non-radial null geodesies from the 
naked central singularity vanishes in the scenario considered here. 
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1 Introduction 



The continual gravitational collapse of a massive matter cloud results either in a black hole 
or a naked singularity, depending on the nature of the regular initial data from which the 
collapse develops. An important development is the emerging realization that at least in the 
spherically symmetric case both these final states seem to occur generically[|l]]. 

The theoretical properties and possible observational signatures of a black hole and a 
naked singularity would be quite significantly different from each other. An immediate 
distinction is, in case the collapse ends in a black hole, an event horizon develops well before 
the occurrence of the singularity, and thus the regions of extreme physical conditions (e.g. 
blowing up densities and curvatures) are hidden from the outside world. On the other hand, 
if the collapse develops into a globally naked singularity, then the energy of the regions 
neighbouring the singularity can escape via the available non-spacelike geodesies paths or 
via other non-geodetic, non-spacelike trajectories to a distant observer. 

In such postulated by some authors@, a huge amount of energy could possi- 

bly be released, in principle, during the final stages of collapse from the regions close to the 
singularity. As suggested in these works, perhaps an enormous amount of energy will be gen- 
erated either by some kind of a quantum gravity mechanism, or by means of an astrophysical 
process wherein the region simply turns into a fireball creating shocks into the surrounding 
medium. One way to study the structure of these extreme regions is to examine the complete 
spectrum of all nonspacelike geodesies through which this energy could escape. Even if a 
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fraction of the energy so generated is able to escape to a distant observer, an observational 
signature could be generated. It therefore becomes quite important and interesting to look 
into this possibility in some detail, and to consider possible observable differences in the two 
scenarios from the perspective of a faraway observer. 

We therefore make such an attempt here to study the appearance of the late stages of 
collapse. While our emphasis in the present study will be on examining the null spectrum, 
which is important, we also need to keep in mind the possibility that emissions could as 
well come through other non-spacelike geodesies or non-geodetic paths, which represent high 
energy particles. 

If the collapse ends in a black hole, the neighboring regions of singularity would appear 
dark, hidden within the event horizon. However, the naked singularity scenario could be 
different in principle. Hence we analyze the possible appearance of these high curvature 
regions when the singularity is naked. In this connection, one of the most extensively studied 
model has been that describing the collapse of inhomogeneous dust, where various features 
arising as collapse end states are well-understood. It is known in this case that both black 
holes and (locally and globally) naked singularities, with either asymptotically flat or a 
cosmological background, could develop depending on the nature of the initial data, which 
is specified in terms of the density and velocity profiles for the collapsing shells at the onset 
of the collapse [Q, |], Hence, these provide an ideal situation to study a problem such as 
above and to study the scenario developing in later stages of collapse. 



4 



It has been shown that in the case of collapse ending in a naked singularity, outgoing 
future directed radial null geodesies (RNGs) families terminate in the past at the central 
naked singularity. The observed appearance of a distant object in the electromegnatic spec- 
trum depends on the behaviour of the full spectrum of null geodesies, including radial as 
well as non-radial trajectories. The behaviour and contribution of non-radial null geodesies 
(NRNGs) could be significant towards the appearence of the singularity. Though the study 
of outgoing non-radial null geodesies terminating at the naked singularity is of theoretical 
interest in its own right, it becomes particularly essential in the context of appearance of 
naked singular regions when one realizes that past calculations of the appearance of a col- 
lapsing star (non-singular regions) depended on the behaviour of non-radial null geodesies 
in the late stages of collapse 0. In fact, in the late stages of collapse the main contribution 
comes from non-radial null geodesies, which form a ring at R = 3M. This ring then de- 
cays exponentially, giving a distinct rapid decline of luminosity within a short time interval. 
It is also of interest to check if timelike non-radial geodesies (NRGs) too came out of the 
singularity. 

We thus need to understand and analyze the behaviour of non-radial geodesies in the 
vicinity of a naked singularity towards such a purpose, which has not been done so far. 
Such an analysis is not only of theoretical interest in its own right, but would also clarify 
several issues such as above, thus revealing the structure of the naked singularity better. 
In the following, Section 2 describes basic equations governing the Tolman-Bondi-Lamaitre 
(TBL) models, and in Section 3 we consider non-radial geodesies emerging from the naked 
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singularity forming in such a collapse. Section 4 discusses the luminosity aspect and in 
Section 5 some conclusions are proposed. 



2 The TBL model and gravitational collapse 

The collapsing dust cloud is described by the TBL metric 0, which is given in co moving 
coordinates as, 



ds z = _ rft 2 + -^jdr 2 + R 2 dQ 2 . (1) 



The energy-momentum tensor is that of dust, 



T ij = edidi, e = e(t,r) = -^, (2) 



where e is the energy density, and the area radius R = R(t, r) is given by 



r 2 = m + (3) 



Here the dot and prime denote partial derivatives with respect to the coordinates t and r 
respectively, and for the case of collapse we have R < 0. The functions F and / are called 
the mass and energy functions respectively, and they are related to the initial mass profile 
and velocity distribution of shells in the cloud. 
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At this point we limit ourself to marginally bound case because of simplicity and clarity 
of the analysis. The results in general case would however be the same. For a marginally 
bound cloud (/ = 0), the integration of equation (3) gives 



where £o( r ) is a constant of integration. Using the coordinate freedom for rescaling the radial 
coordinate r we set, 



At the time t = t (r), the shell labelled by the coordinate radius r becomes singular where 
the area radius R for the shell becomes zero. We consider only the situations where there 
are no shell-crossings in the spacetime. A sufficient condition for this is the density be a 
decreasing function of r, which may be considered actually to be a physically reasonable 
requirement, because for any realistic density profile the density should be higher at the 
center, decreasing away from the center. The ranges of coordinates are given by, 




(4) 



R(0,r)=r, 



(5) 



which gives, 




(6) 



< r b < oo, — oo < t < t (r), 



(7) 
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where r = r b denotes the boundary of the cloud. The quantity R', which is also needed later 
in the equation of the geodesies to check the visibility or otherwise of the central singularity, 
can be written as, 



R' 



F'R 



rF' i 

+ (i-^)v^- 



(8) 



3F 



3 Non-radial geodesies in TBL models 



The regular center in a spherically symmetric spacetime can be the source of only radial 
geodesies. So when the central singularity is covered, the center is observed by a distant 
observer by intercepting radial null geodesies till the time the event horizon forms, beyond 
which nothing is seen. An outside observer would see radial photons from the center, which 
would disappear much before the formation of the singularity. These photons do not stay 
in the cloud longer. At this point it is important to note that in earlier discussions of a 
collapsing star entering the event horizon, the role of non-radial geodesies cannot be over 
emphasized. The non-radial null geodesies emitted in late stages of collapse tend to revolve 
around the star and thus take longer time to reach the distant observer. Therefore, even after 
the star has gone through it's Schwarzschild radius, it leaves a ring of photons at R = 3M. 
So a distant observer observing late stages of collapse sees a bright ring at R = 3M which 
decays exponentially. In the context of formation of naked singularity, therefore, it is quite 
essential that we examine this issue carefully and investigate the non-radial spectrum of null 
geodesies. 
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The tangents to the outgoing geodesies are given by 



r w dt P r _ rfr VP 2 - 1 2 + 5i? 2 
dA R' d\ RR> 



(9) 



fT = 



def) I 
d\ = R 2 ' 



(10) 



Where P satisfies the differential equation, 



R R' 



R 2 RR> 



- VP 2 -l 2 + BR 2 — + BR = (11) 
R 2 



We can work in the equatorial plane 9 — | , with K e = (dO/dX) = 0, and due to spherical 
symmetry, this could then be rotated to recover the same qualitative features for the full 
spectrum of null geodesies. In the above equations B = for null geodesies, B = — 1 for 
timelike geodesies, and B = 1 for spacelike geodesies. 

These equations can be written in the (r, R) plane as, 



dR R' 



dr a ar a 1 



1 - 



R y/l ~ l 2 /P 2 + BR/P 2 



(12) 



where a is a constant fixed by demanding that R'/r a 1 remains finite. For the quantities <fi 
and P we have, 



I 



dr P Ry/P 2 -l 2 + BR 2 ' 



(13) 
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and, 



dP 

dr 



PR' 



VP 2 -l 2 + BR 2 



1 [F\3R' 



F' 



+ 



BRR' [F 



(14) 



R 




F 



y/P 2 - P + BR? V R 



respectively. 

From the above equations it follows that for t < io(0) and r = 0, the tangent equations 
give for the self-consistency I = 0. Thus no NRGs would be radiated by the center before 
the singularity is formed at the center at r = 0. 

We now turn our attention to non-radial geodesies in TBL models. Defining u = r a , 
X = R/u, and using the l'Hospital's rule we can write, 



where U(X, u) = dR/du (along the geodesies). If the above equation has a real positive root 
X = Xo then the singularity is at least locally naked ||. We want to solve these equations 
in a self-consistent manner near the center, requiring that the outcoming null geodesies have 
a well-defined tangent at the singularity in a suitable plane. For simplicity and clarity, we 
assume the mass function to have the form, 



X = lim X 




lim — = U(X , 0) 

^0,u^0 du 



or, 



U(X,0)-X = V(X) = 0. 



(15) 
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F = F r 3 + F n r 3+n + higher order terms, 

and we first consider only the null geodesies (B = 0). Later we generalize the results to 
other geodesies. 

Now, we can see in n < 3 cases, assuming that l/P < 1 and checking for self-consistency, 
that there are two kinds of behaviors possible for R, as given by, R ~ X r 1+2n ^ 3 (where 
X = (--F„/2P ) 2/3 ) or R « X r 3 « P r 3 || |]. In the first case, from equation (0) 
(assuming that P remains nonzero) we get, 

P w p Q r a = P r 1+2n/3 , (16) 

where Pq is constant of integration. So we get a contradiction, because 1 — I 2 / P 2 — > — oo. It 
follows that we cannot have NRNGs coming out of the singularity along this direction. 



In the second case (R ~ X r 3 w F o r 3 )0, ^] assuming, 1 — I 2 / P 2 > we get for equation 
(0), 



P = P e 6 ( 3 -) F o /2 . (17) 



Then, in the limit to the singularity P blows up exponentially so we get 1 — I 2 / P 2 — > 1 and 
we have a self-consistent solution for the non-radial geodesies. One can also check that there 
are no self-consistent solutions for NRNGs along any other directions. 
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In the case n = 3, we can uniquely fix a = 3 and assuming, along the geodesies R = X r 
we get from (|T4D, 



3(2tg-2h -l) 



P « P r 26 o(6o-D , (18) 



where b = JX /F . For P to blow up we need (26q — 2b — 1) > 0. This gives, 



1 + V3 



b = VXo/Fo < b 0crit = (19) 

Assuming that P blows up, the root equation ( |T5D reduces to the same equation as for the 
radial null geodesies case. Introducing X = F x 2 and £ = F 3 /F^ 2 the root equation then 
becomes, 

2x 4 + x 3 - £x + f = 0. (20) 

From the theory of quartic equations, this admits a real positive root for £ < £ cri4 = —(26 + 
15/v / 3)/2. Basically, if this condition is satisfied then we have two real positive roots. The 
smallest value for the larger root x\ is the same as the largest value of the smaller root xi and 
it is same as b 0crit , and this is achieved at £ = ^ crit . So the condition ( |19"D is never satisfied 
along the larger root x\ (Cauchy horizon), and is always satisfied along the smaller root x^- 
So P — > along the larger root, and we cannot get self-consistent real solutions for the set of 
equations for the null geodesies. That means we cannot have NRNGs coming out along the 
Cauchy horizon direction. Along the smaller root direction this condition is always satisfied, 
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P blows up, and we have a self-consistent solution for the geodesic differential equation. 
That means we can have NRNGs coming out of the singularity only along the direction of 
the smaller root. 



If n > 3 we have to fix a = 1 + 2n/3 > 3, so F/R — > oo, and we cannot have any null 
geodesies coming out of the singularity. 

Now we will check what happens for timelike (B = —1) and spacelike (B = 1) NRGs. 
Even in these cases, in equations ( plfflf) , for all the singular geodesies BR 2 goes to zero as 



we approach the center (r = 0) and it is always negligible compared to the other terms. So 
provided P remains nonzero near the center, the root equation (|T"5|) remains the same as that 
for the null geodesies. Now if we also solve the equation for P (O) in self-consistent manner 



we see that for the smaller root P (and so also K l ) blows up and we get a self-consistent 
solution for the geodesic equations near the singularity, and the behavior of these geodesies 
near the singularity is very similar to the null geodesies. Similarly we can easily check that 
no non-radial geodesies can come out along the larger root direction. 

We now consider the global behavior of non-radial null geodesies coming out of the 
singularity in various cases, and check whether there are any NRNGs which can be actually 
seen by a faraway observer. Basically, we try to check for singular NRNGs (outgoing non- 
radial geodesies which terminate in the past at the singularity), and examine with what 
maximum value of I (I — l ma x) can they meet the boundary of the cloud r = r c , when the 
boundary has area radius R = R c . Here we fix the boundary in such a way that the density 
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goes to zero smoothly at the boundary. The outside solution is Schwarzschild, with total 
mass m = F(r c )/2. So with our condition, we fix, 



Figures la, 2a, show the graphs of R c versus l m ax/P c - (Where P = P c at R = R c and 
we give this as the initial condition.) We see that l max / 'P c increases as we decrease R c . This 
is in a way expected as the geodesies with larger value of l ma x/P c revolve more around the 
center, i.e. they remain in the inner region for longer times, so they are more likely to get 
trapped before reaching the larger value of R. Figures lb and 2b, show <p verses r while 
Figures lc and 2c show R verses r for some singular and non-singular geodesies. 

Considerations above point to NRNGs terminating in the past at a naked singularity 
and could reach a distant external observer. Unlike the regular center of the cloud, which 
cannot have non-radial geodesies terminating there, we have a distinctly new scenario where 
an observer could intercept NRNGs from the naked central singularity. Because of the high 
curvature in the singular regions non-radial photons revolve around the center and stay for 
a longer time in the cloud before reaching the surface and hence to the distant observer. 
Therefore unlike the black hole case, where the center of the cloud disappears rather quietly, 
in the case of a naked singularity, the center would appear as nothing unusual upto the 
time singularity is formed. However, once the singularity is formed, it would appear rather 
differently. 

Let us consider this phenomena in some detail. For a distant observer, till the receipt 




3 +nF r 
~^~F 



n\l/n 
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0.12 




Figure 1: Various plots for F = r 3 — 4r 5 case, (a) Plot showing the maximum value of 
Imax/Pc for various values of R c . (b) Plot showing the geodesies with I = l max in (r, 0) plane 
for two different values of R c . (c) Same geodesic in (r, R) plane. 
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Figure 2: Various plots for F = r 3 — 30r 6 case, (a) Plot showing the maximum value of 
Imax/Pc for various values of R c . (b) Plot showing the geodesies with I = l max in (r, 0) plane 
for two different values of R c . (c) Same geodesic in (r, R) plane. 
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of the first singular RNG from the naked singularity, the regular center would be observed 
through the RNG and would appear as a regular center inside the cloud. After the formation 
of naked singularity, however, the observer would start also intercepting NRNGs which arrive 
later than their RNG counterpart. The observed position through the NRNGs depends on 
the value of the impact parameter I. Higher the value of the impact parameter more time 
is spent by these photons in revolving inside the cloud and therefore taking longer time 
to come out. The maximum value of the impact parameter / max (t ) is a function of time. 
Therefore, at any instant of time, after receiving the first RNG from the naked singularity, 
the observer would receive NRNGs within the range l max (t) > I > (for the cases when 
there are families of singular RNGs coming out). To a distant observer, the photons would 
appear to be coming not only from the center but from a spherical region (with r = as 
center) having a radius equal to l max . In this sense therefore, the center would appear as a 
spherical ball with an expanding radius. The expansion continues till a maximum value of I 
is reached. Apart from the visible appearance of the naked singularity in the electromagnetic 
spectrum, the observations in terms of high energy particles traveling along non-spacelike 
curves would also be the same. That is, unlike the black hole scenario, in case of a naked 
singularity forming, the timelike particles would also appear to be coming from an expanding 
ball. 
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4 Luminosity function for the singularity 



Normally most of the luminosity (energy) of an object comes from the non-radial rays. We 
have shown above that an infinite number of non-radial geodesies come out of the naked 
singularity, along the non-Cauchy horizon direction. Only one RNG comes out of the sin- 
gularity along the Cauchy horizon direction (larger root) p. So, if any radiation comes out 
of the singularity, it will be normally expected to come out along the smaller root direction 
through the rays with various values of the impact parameter I. From such a perspective, 
we need to consider the luminosity function for these rays for various observers. 

Let us consider an observer at comoving coordinate r = r . The observed intensity I p of 
a point source is 

'^aJtW (21) 

Where P is the power radiated by the source into the solid angle SQ, and A is the area 
sustained by the rays at the observer. The redshift factor (1 + z) 2 appears because the power 
radiated is not the same as power hitting the area at observer. In case of RNGs Aq oc Rq 
where Rq is the radius of the observer. 

Let -u a ( s ) and w a ( ) be the four-velocities of the source and the observer and let E\ and E 2 
be two events connecting the source and the observer through the RNG. The redshift factor 
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is given by 



l + z = 



[KgU a {s) } El 



(22) 



where the numerator and denominator are evaluated at events E\ and E 2 , at the source and 
observer respectively, with 



Taking the source as the naked center at r = and the observer at r = r we have 



In the evaluation of the redshift, the behavior of the tangent vector component K* is 
important. It is finite at the nonsingular observer who is sufficiently far away from the 
center, i.e. r Q » 0. Therefore, the behavior of K l at the naked singularity determines 
basically the behavior of the redshift factor. As discussed earlier, along the trajectories of 
interest K l diverges very rapidly (exponentially if a < 3, and by a large power law if a = 3). 
This means that the redshift diverges very rapidly for the singular rays of our interest and 
for any non-singular observer it will be infinite. (In a way the exponential divergence of the 
redshift for a < 3 is expected as these rays stay very close to the apparent horizon near 
the singularity and in such a situation the redshift in the Schwarzschild case also diverges 
exponentially.) As the redshift diverges, and A is finite, classically the luminosity of the 
naked singularity along such families should vanish. 




u\o) = St- 



(23) 
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It is very difficult to get rid of this dominance of rapidly diverging K l near the singularity 
on the redshift z, and so also the dominance of this diverging redshift on the luminosity 
That means the luminosity of the naked singularity along such NRNGs for an observer will 
be zero, i.e. no energy will reach the observer along such families from the naked singularity, 
at least classically. Therefore, in the case above, the naked singularity may not be physically 
visible to faraway observers directly by means of emitting light. However, for the a < 3 
naked singularity cases, if an enormous amount of radiation is emitted just along the Cauchy 
horizon, it may be possible that a measurable fraction of the emitted energy will reach an 
observer. This may happen as along this wavefront K f , and so also the redshift remains 
finite 0. But this observation of luminosity would only be instantaneous as only a single 
light ray is allowed to escape. 

All the same, the possibility of mass emission via timelike or non-spacelike non-geodetic 
families of paths coming out from the naked singularity remains open. In the case of such 
a violent event being visible, particles escaping with ultrarelativistic velocities cannot be 
ruled out from this neighbourhood. It is also to be noted that the classical possibilities 
such as above regarding the probable light or particle emission, or otherwise, from a naked 
singularity may not perhaps offer a serious physical alternative one way or the other. The 
reason is, in all physical situations, the classical general relativity would break down once 
the densities and curvatures are sufficiently high so that quantum or quantum gravity effect 
become important in the process of an endless collapse. Such quantum effects would come 
into play much before the actual formation of the classical naked singularity, which may 
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possibly be smeared out by quantum gravity. The key point then is the possible visibility or 
otherwise, of these extreme strong gravity regions, which develop in any case, in the vicinity 
of the classical naked singularity. It is then the causal structure, that is, the communicability 
or otherwise, of these extreme strong gravity regions that would make the essential difference 
as far as the physical consequences of naked singularity formation are concerned. 

In the black hole case, resulting from the collapse of a finite sized object such as a massive 
star, such strong gravity regions, or what we may call 'fireballs', will be necessarily covered 
by an event horizon of gravity, well before the curvature conditions became extreme (e.g. 
well before the collapsing cloud went to the Planck size). In such a situation, the quantum 
effects, even if they were causing qualitative changes closer to the singularity, will be of no 
physical consequences, because no causal communications are allowed from such regions. 
On the other hand, if the causal structure were that of a naked singularity, communications 
from such a quantum gravity dominated extreme curvature ball would be visible in principle, 
either directly, or via secondary effects such as shocks produced in the surrounding medium. 
Then we may have a chance to observe directly the quantum gravity effects from such fireballs 
generated due to stellar collapse. 

5 Conclusions 

Studying NRNGs we have demonstrated that whenever the singularity is naked, along with 
the RNGs, NRNGs also come out of the singularity. One can also say that if NRNGs came 
out of the singularity then RNGs will also come out, because we have shown that eventually 
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the root equation for the existence of geodesies does not depend on the value of the impact 



parameter I. This is similar to the recent paper ]I0| on null geodesies, which came while 
this was being written. We have also shown that similar results exist for the timelike and 
spacelike geodesies also. Basically we see that in comoving coordinates all the geodesies have 
a very similar behavior near the central singularity. 

The existence of non-radial geodesies coming out can in a way make the naked singularity 
look like having a finite area to an outside observer. This possibly happens because even 
though the singularity is at the center of symmetry and geodesies can come out of it, gravity 
is very powerful and dominant for such trajectories. 

We have also shown that NRNGs come out only along the direction of the smaller root, 
i.e. no NRNG comes out along the Cauchy horizon direction which corresponds to the larger 
root. Together with the earlier results H, this means that only a single radial null geodesic 
comes out of the naked singularity along the Cauchy horizon direction. Studying the global 
behavior of geodesic families using numerical methods, we have shown that even when near 
the singularity any value of / is allowed, only the geodesies with a certain maximum value 
of lmax/Pc can reach any given outside observer. This value depends on the position of the 
observer and is larger if the observer's area radius is smaller. This is expected as geodesies 
with larger value of / will stay near the center for a longer time, and as the cloud is collapsing 
they are more likely to get trapped. The numerical study also shows that typically all such 
geodesies undergo finite number of revolutions while they go out. The main reason for this 
is that effects of gravity are very dominant near the center. Typically these geodesies can 
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revolve around the center at the most a few times before escaping. 

Studying the redshift and luminosity along the various geodesies we have shown that, 
apart from one special RNG (the Cauchy horizon) in a < 3 case, along all other singular 
trajectories the redshift diverges for any comoving (non-singular) observer and so the lu- 
minosity reaching from the naked singularity will be zero along such families. One could 
argue that even if a single photon, or a single wave front carrying huge energy escapes from 
the singularity with a finite redshift, that may destroy the cosmic censorship, because that 
can alter the qualitative picture considerably. However, normally we do not expect a single 
wavefront to emit an arbitrarily large amount of energy, though one does not know what 
happens near such extreme regions. What we may say then is, for a = 3 (i.e. n — 3) naked 
singularity energy along null geodesies is censored for all the trajectories. For a < 3 (i.e. 
n < 3) it is censored apart from the first trajectory, i.e. for trajectories along the larger 
root direction. Thus for regions of spacetime where the curvature diverges fast enough, the 
redshift is infinite thus censoring energy. In a sense this result can be considered to be 
supporting the cosmic censorship hypothesis, if we mean by the later a statement such as 
general relativity allows the occurrence of naked singularities, however, they may not directly 
radiate away energy to outside observers. What this means exactly is that though general 
relativity allows the occurrence of naked singularity, the radiated energy in the electromag- 
netic spectrum does not reach the distant observer, at least in the dust case. We need to 
check such a statement for timelike and non-spacelike non-geodetic paths coming from the 
naked singularity, and also for equations of state other than dust, when the naked singularity 
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may have a complicated topology. We note of course, that cosmic censorship is not really 
a statement about the energy escape, it is essentially a basic postulate about not having 
outgoing causal curves from the singularity. 

Though for simplicity and clarity we have shown these results for marginally bound case, 
it would be possible to generalize the same to non-marginally bound cases using a similar 
method, and depending on the value of a needed to make R! /(r^ -1 ) finite, the results will 
be similar. This will be the case as depending on the value of a various functions involved 
in this analysis can be expanded in the similar way near the central singularity. Further, we 
can expect the timelike NRGs as well to come out of the naked singularity even in these cases 
as the study shows that they also have very similar behaviour to that of the null geodesies 
in the vicinity of the singularity. 
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